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Abstract 

To study orthogonal arrays and signed orthogonal arrays, Ray-Chaudhuri 
and Singhi (1988 and 1994) considered some module spaces. Here, using a 
linear algebraic approach we define an inclusion matrix and find its rank. In 
the special case of Latin squares we show that there is a straightforward al- 
gorithm for generating a basis for this matrix using the so-called intercalates. 
We also extend this last idea. 

Keywords: Orthogonal arrays, Latin squares, basis for inclusion matrix, Latin 
trades 



1 Introduction and preliminaries 

To show the existence of signed orthogonal arrays, Ray-Chaudhuri and Singhi 
considered a space of linear forms in variables and calculated its rank, see ff3l . 
Later they pointed out an error in their calculation and provided a correction, 
see fl4l . Here we define a natural inclusion matrix corresponding to orthogo- 
nal arrays and signed orthogonal arrays. We compute the rank of this matrix and 
study bases of its null space. This provides helpful insight for studying these ob- 
jects. In the special case of Latin squares we show that there is a straightforward 
algorithm for generating a basis for this matrix using the so-called intercalates. 
We also extend this idea for more general cases. 

We follow the notations of 1 13 1 as much as possible. Let V := {0, 1, . . . ,v — 
1} and V be the set of all ordered fc-tuples of the elements of V, i.e., V k :— 

{(xi,...,Xk) | Xi € V,i = l,...,k}. Also, let V} := {(ui, . . . , Ut)i \ u t € 
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V.i = 1, . . . , t}, where I is a subset of size t of the set {1, . . . , k}. For a pair of 
elements of V k and Vf, where / = {ii, . . . ,i t } and i\ < ■ ■ ■ < i t , we define: 
(ui, . . . ,Uf)j e (xi, . . . ,x k ) uj=x 4j , j = l,...,t. 

An orthogonal array OA t (u, fe, A) on a set F is a collection of fc-tuples of el- 
ements of V such that for each / C {1, . . . , k}, \I\ = t, every element of Vj 
belongs to exactly A elements of the collection. Orthogonal arrays were first de- 
fined by Rao |[T2l and have been used in studying designs and codes. 

Next we define the ^-inclusion matrix M(t-(v, k)). Columns of this matrix cor- 
respond to the elements of V k (in lexicographic order) and its rows correspond to 
the elements of UjVj, where the union is over all t-subsets of {1, ... , k}. The en- 
tries of the matrix are or 1, and are defined as follows. Mr Ul! „^ u A I / Xli „^ Xk ) = 

1 (ux,...,U t )l S (xx,...,Xk). 

Example 1 In Figure\I\ the matrix M(2-(3, 3)) is shown. To make it more read- 
able, entries of the matrix are represented by "■ " signs. The "*" signs may be 
ignored for time being, they will be explained in Section\2\ 
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Figure 1: M(2-(3,3)) 
A Latin square L of order v is a v X v array with entries chosen from a set, 
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say V — {0, 1, . . . , v — 1} in such a way that each element of V occurs pre- 
cisely once in each row and in each column of the array. For ease of exposi- 
tion, a Latin square L will be represented by a set of ordered triples Lij) 

element occurs in cell of the array}. 

It is folkloric that any Latin square of order v is equivalent to an OA2(v, 3, 1). It 
is easy to see that any OAt (v, k, A) can be thought of a solution to the equation 



where M = M(t-(v,k)), 1 is a vector of appropriate size with all components 
equal to 1, and F is a non-negative integer-valued frequency vector, i.e., F(x) 
represents the number of times that OA contains the ordered fc-tuple x. Our dis- 
cussion will be based on the field of real numbers. 

In Section [2] we find the nullity and the rank of M(t-(v, k)). In Section [3] we 
show that a very simple basis exists in the case of M(2-(u, 3)), i.e., when the OA 
corresponds to a Latin square, which consist of so-called intercalates. In Section|4] 
we generalize the result of Section[3]for any M(t-(v, t + 1)). 

2 Orthogonal arrays 

The main result of this section is the following theorem, 
Theorem 1 The rank of the matrix M.(t-(v, k)) is equal to 



This theorem results from the following lemmas. But first we need the following 
notations. For every ordered fc-tuple x = {x\ ,...,£&), the set F x is defined as 



Also, we define A x = {i \ Xi ^ 0}, L x = \A X \, and let C x denote the column of 
the matrix M corresponding to the fc-tuple x. 

Lemma 1 For every x G F y , x =/= y, we have L x < L y and x -< y, where -< 
denotes the lexicographic order. 

Proof. Clearly, for every i, if Xi / 0, then j/j = X{ ^ 0. Therefore, L x < L y . 
Now if the equality L x = L y holds, then for every non-zero yi, the corresponding 
Xi is non-zero and hence equal to j/j. This implies x = y, contradicting the 
hypothesis. For the second part, notice that for every i, either Xi — yi, or = 
Xi < yi. Hence x^y and since x ^ y, we have x -< y. ■ 



MF = Al, 



(1) 




i=0 



-Fx = {(zi, ■ ■ ■ , Zk) | Zi G {0, Xi}, i=l,...,k}. 
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Lemma 2 The number of linearly independent rows ofWL(t-(v, k)) is at least the 
number of columns C x with L x < t. 

Proof. We show that for every column C x with L x < t, there exists a correspond- 
ing row such that its pivot 1 (that is the first 1 in that row) is in column C x . These 
rows are clearly linearly independent. Corresponding to the column C x , we con- 
struct an element u/ = (ui , . . . , utji of V} as follows. Since L x < t, we can pick 
a set / = {ii, . . . , it} with A x C I. Note that there may be more than one choice 
for I, but choosing any of them will serve our purpose. For every j = 1, . . . , f, 
we let Uj = Xj,. This defines an element Uj such that Uj 6 x. Therefore, there 
is a 1 in the intersection of the column C x , and the row corresponding to Uj in 



Next, we prove that every y such that u/ G y satisfies x ^< y. This would 
imply that the first 1 in the row corresponding to u/ lies in the column C x , which 
completes the proof. Since for every j = 1, . . . , t, Uj = x^ and Uj = y^ , we 
have Xi = yi for every i g I. Also, since Xi = for every i g I, we have 
Xi G {0,yi} for every i, or in other words, x G F y . Therefore, by Lemma [T] 
xr<y- ■ 

Lemma 3 rank(M(t-(«, fc))) > ^Lo Cl)( v - 1 Y- 

Proof. By Lemma 12 the rank of M is at least the number of columns C x with 
L x < t. It is easy to note that the number of vectors x for which L x — i is equal 
to ( .) (v — iy. This completes the proof of the lemma. ■ 

As it can be noted in Figure [U in each column C x where at least one of the 
components of x is 0, i.e. L x < 2, there exists at least one row having a pivot 1 
in that column. For example for the column Coio both (0, l){i,2} an d (lj 0){2,3} 
rows have such property. All such columns and one of the rows corresponding to 
that column are indicated with a "★" sign and such 1 's are shown in bold face and 
underlined. 

To show the other direction in Theorem[U we prove the following lemma. 
Lemma 4 For every vector x G V k with L x > t, we have 



Where is a vector of appropriate size with all components equal to 0. 

Proof. It is enough to focus on a fixed row, say Uj, and count the number of ones 
in the intersection of this row and columns C y for y G F x , by taking the signs in 



M(t-(v,k)). 




(2) 
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the above expression into account, and show that the corresponding entry in the 
left-hand side of the above equation is 0. 

Consider a row uj with / = {ix, . . . ,i t } where i\ < • • • < i%. 

If there is an element ij E I \ A x such that Uj ^ 0, then all the entries in the 
intersection of this row and columns C y with y E F x are 0, since for all such y, 
Ui. is and therefore is not equal to Uj. Thus, the entry in the left-hand side of 
Equation (O in the row corresponding to uj is 0. 

Now, suppose Uj — 0, for every j where ij E I \ A x . Consider the set Y of 
vectors, y E F x , such that the entry in the intersection of the row corresponding 
to uj and the column C y is 1. For every y € Y, a = \A X n I\ of its entries have a 
value equal to the corresponding entry in u/. Therefore, there are i x — a entries 
in y that can take either a value of 0, or the value of the corresponding entry in x. 
We call these L x — a entries the free entries of y. Consider the set of all y 6 Y 
that have j non-zero free entries (i.e., are equal to the corresponding value in x). 
The number of such y's is and for each such y, L y — j + a — £, where £ 

is the number of zeros in the intersection of / and A x . Therefore, the entry in the 
row corresponding to Uj in the left-hand side of Equation (f2j) is equal to 

L gVl)^-< ( L ~- a ) = (_i)-C "ff (_!)* ( L * - a \ = o. 

j=Q \ J J - =Q \ J / 

Thus, all entries of the vector in the left-hand side of Equation (O are 0. ■ 
Lemma 5 rank(M(t-(u, k))) < J2l =0 ( ■) [v ~ l) 1 - 

Proof. By applying Lemma|H repeatedly, to any column C x with L x > t, we can 
write C x in terms of columns C y with L y < t. Thus the latter columns form a 
spanning set for the column space of M(t-(v, k)). We noted earlier that there are 
exactly ^j_ (f) (v — l) 1 of such columns. ■ 

Theorem[T]follows from Lemma|5]and Lemma[3] 

3 Latin squares and Latin trades 

As we noted earlier, a Latin square of order v may be viewed as an OA2(v, 3, 1). 
So the matrix M(2-(i>, 3)) is of special interest. In this section we find a basis for 
M(2-(w, 3)) which consists of the so-called intercalates. 

We start with a few definitions. A partial Latin square P of order v is a vxv array 
in which some of the entries are filled with elements from asety = {0,l,...,u— 
1} in such a way that each element of V occurs at most once in each row and at 
most once in each column of the array. In other words, there are cells in the 
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array that may be empty, but the positions that are filled conform with the Latin 
property of array. Once again a partial Latin square may be represented as a 
set of ordered triples. However in this case we will include triples of the form 
0) and read this to mean that cell of the partial Latin square is empty. 
The set of cells Sp = | (i,j;Pij) E P, for some Pij G V} is said to 

determine the shape of P and \Sp\ is said to be the volume of the partial Latin 
square. That is, the volume is the number of nonempty cells. For each row r, 
< r < v — 1, we let lZ r p denote the set of entries occurring in row r of P. 
Formally, lZ r p = {P r j \ P r j 6 V A (r, j; P r j) e P}. Similarly, for each column 
c, < c < v - 1, we define Cp = {P lc \ P lc e V A (i, c; P lc ) £ P}. 

A Latin trade, T = (P, Q), of volume s is an ordered set of two partial Latin 
squares, of order v, such that 

1. S P = S Q: 

2. for each G S P , P^ / Qij, 

3. for each r, < r < v — 1, TV P = TZ-q, and 

4. for each c, < c < v - 1, Cp = C C Q . 

Thus a Latin trade is a pair of disjoint partial Latin squares of the same shape and 
order, which are row-wise and column-wise mutually balanced. We refer to the 
shape of a Latin trade T as the shape of the individual components P and Q. 

Example 2 Below is an example of two partial Latin squares which together form 
a Latin trade of order 5 and of volume 19. To conserve space we will display a 
Latin trade by superimposing one partial Latin square on top of the other, and 
using subscripts to differentiate the entries of the second from those of the first, as 
shown below. 
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Figure 2: A Latin trade 

The concept of a Latin trade in a Latin square is similar to the concept of a mu- 
tually balanced set or a trade in a block design, see 0. The same as trades in 
design theory, the discussion of Latin trades is related to intersection problems. 
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For example, they are relevant to the problem of finding the possible number of 
intersections for Latin squares (see (7), 0, E), and HI). Also Latin trades arise 
naturally in the discussion of critical sets in Latin squares (see for example IfTOl 
and O]). 

Latin trades have been studied by many authors. Fu and Fu [6| used the term 
"disjoint and mutually balanced" (DMB) partial Latin squares, Keedwell [ 10 1 used 
"critical partial Latin square" (CPLS), while Donovan et al. J4| used the term 
"Latin interchange". Adams et al. l suggest the terminology '2- way Latin trade" 
for consistency with similar concepts in other combinatorial structures such as 
block designs, graph colouring, cycle systems, etc. See for instance (9), 021, El, 
and for further use of trades. 

A Latin trade of volume 4 which is unique (up to isomorphism), is called an 
intercalate (see Figure©. 
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Figure 3: An intercalate 

Similar to orthogonal arrays which correspond to solutions of the Equation (fTJ, it 
is clear that any Latin trade may also be treated as a solution T, to the equation: 

MT = 0, (3) 

where M — M(2-(i>, 3)) and T is a (signed) frequency vector derived from the 
trade T = (P, Q), i.e., 

fl if (i,j;k)eP 

T ijk = l-l if (i,j;k)eQ 

[0 otherwise. 
Therefore Latin trades are in the null space of M(2-(w, 3)). 

Theorem 2 There exists a basis for the null space of the matrix M(2-(i>, 3)) con- 
sisting only of intercalates. 

Proof. Latin trades are in the null space of M(2-(v, 3)). By Theorem[T]we know 
that 

null(M(2-(u,3))) = {v-lf. 
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For each k; 1 < k < v — 1, consider the ij'fc'th intercalate defined in 
Figure 4. 
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Figure 4: Basis intercalates 

There are (v — l) 3 of them. The vectors corresponding to these intercalates are 
independent, as for example, the frequency vector of the ijk'fh intercalate has 
an entry —1 in the k) coordinate while all others have in that coordinate. 
Therefore, they form a basis for the null space of M(2-(v, 3)). ■ 

The above theorem shows, existentially, that every Latin trade can be written as the 
sum of intercalates. In [5|, Donovan and Mahmoodian introduced a simple com- 
binatorial algorithm which enables one to compute such a decomposition. But by 
linear algebraic approach and knowing a basis which consists only of intercalates, 
makes it straightforward to do this task. We give an example of this method: 
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Figure 5: An example of the algorithm 



4 A basis for the null space of M(t-(v, t + 1)) 

In this section we generalize the notion of Latin trades and find a basis for the null 
space of M(t-(v, t + 1)). First we note that each Latin trade T = (P, Q) may 
be represented by a homogeneous polynomial of order 3 as follows. The polyno- 
mial is over a non-commutative ring, hence the terms are ordered multiplicatively 
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(meaning that x^x^x^ is different from, say x^x^x^): 

P^Xq , X\ , . . . , X v — \ ) — ^ ^ %iiX{ 2 Xi 3 ^ ^ ^31^32^33' 

Note that the positive terms correspond to the elements of P while the negative 
terms correspond to the elements of Q. For example the intercalates which form 
a basis for the null space of M(2-(i>, 3)) in Theorem[2]are: 

P(x ,Xi, . . . ,X v -i) = X X X + X XjX k + XiX X k + XiX 3 X 

- x x Xk - xqXjXq - XiX xo - XiXjXf. 
= (x - Xi)(x - Xj)(x - x k ), l<i,j,k<v-l. 

Similarly, each homogeneous polynomial P(xq, x%, . . . , x v ~i) of order fc, whose 
terms are ordered multiplicatively, corresponds to a frequency vector T. If T 
satisfies 

MT = 0, (4) 

where M = M(t-(v, k)), we call it a t-(v, k) Latin trade. So, a Latin trade de- 
fined in Section [3] is also 2-(v, 3) Latin trade. Any t-(v, t + 1) Latin trade of the 
following form will be called a t-(v, t + 1) intercalate: 

P(X(), X\ , . . . , X v — i) — {^ii ~ x ji ) (^42 x 32 ) ' ' v^it+l X 3t+l )i 

where i m and j n G {0, . .. ,v — 1}, and for each /, ii is distinct from ji. A 
2-(v, 3) intercalate is simply an intercalate in a Latin square, defined in the previ- 
ous section. 

Theorem 3 There exists a basis for the null space of the matrix M(t-(v, t + 1)) 
consisting only oft-(v,t + l) intercalates. 

Proof. It can easily be checked that t-(v, t+ 1) intercalates are included in the null 
space of M(t-(v,t + 1)). By Theorem[T]we know that null(M(t-(u, t + 1))) = 
(v-l) t+1 . 

Consider the following set of (v — 1)* +1 intercalates: 

P(x ,xi, . . . ,x v -i) = (x - x il )(x - x i2 ) ■■■(x - x it+1 ), 

1 < ii,i2,---,ik < v- 1- 

The intercalates in this set are independent. For example, the frequency vector 
of the iii2 • ■ -it+i'th intercalate has a non-zero entry, namely (— 1)* +1 , in the 
(ii,i2,. ■■ , «t+i)th coordinate, while all others have in that coordinate. There- 
fore, they form a basis for the null space of M(t-(v, t + 1)). ■ 

Finally, we note that given the above generalization of the concept of Latin trades, 
many questions similar to the ones in the theory of i-trades in block designs, may 
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be raised. For example, it would be interesting to characterize the possible support 
sizes of t-(v, k) Latin trades. 
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